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ANGLE-OF-ATTACK ANALYSIS FOR PROJECT FIRE 1 
PAYLOAD REENTRY FLIGHT 


By Gerard E. Woodbury 
Langley Research Center 

SUMMARY ^ / O 3^ 

An investigation of the flight motions of the Project Fire 1 payload was made through 
the use of digital- computer simulation techniques since limited motion instrument data 
did not permit more direct methods of analysis. The simulations were obtained from a 
six-degree- of-freedom computer program which numerically solved the equations of 
motion. The complete nonlinear wind-tunnel data and calculated time histories of the 
mass characteristics of the body were used in the program. The only data available for 
simulation consisted of continuous data from the yaw-rate gyro and commutated acceler- 
ometer data. Some of the rate- gyro data were lost, however, since the excursions of the 
rate gyro oscillations exceeded calibration limits. Also, periodic fadeout of the telemetry 
system caused a loss of both rate- gyro and accelerometer data. A separate study was 
necessary to obtain the initial conditions for the simulations. This study consisted of a 
curve-fit program which used the angular- rate-history solution to the linearized equations 
of motion to approximate the pitch rate from measurements of the yaw- rate gyro. Since 
the resultant force and resultant pitching- moment coefficients of the wind-tunnel aero- 
dynamics were linear up to angles of attack of 35°, this curve-fit program could be used. 
With the initial angular velocities and accelerations obtained at a selected time from the 
curve-fit program, all the remaining initial conditions, with exception of the trajectory 
quantities, necessary to start an exact solution on the six-degree-of-freedom computer 
program were easily determined. The purpose of this investigation was to determine the 
angle-of-attack and motion histories of the reentry package during the heating data period 
of the Project Fire 1 experiment. The reentry package experienced a disturbance at 
about 1666 seconds of elapsed flight time at an altitude of approximately 69 800 meters 
and a velocity of about 11 500 meters per sec. Prior to the disturbance, the body was 
coning with a half-angle of about 1°. 

The results of the analysis indicated that the angle-of-attack envelope increased to 
about 33° after the disturbance. The factors which caused the disturbance could not be 
isolated. The angle-of-attack envelope damped from 33° to about 13° at maximum 
dynamic pressure. The roll rate was initially 18.85 radians /sec. The motion patterns 



of the body were typical of those for bodies which are statically stable and have roll-to- 
pitch inertia ratios greater than 1. The simulations indicated that trim was small. 

INTRODUCTION 

Project Fire is a National Aeronautics and Space Administration research effort to 
obtain heating data on a blunt Apollo shaped reentry body at hyperbolic velocities. The 
first flight test, Fire 1, was launched 14 April 1964, from Cape Kennedy, Florida. 
Approximately halfway through reentry the payload experienced a large unexpected dis- 
turbance which caused severe motions to ensue. Reference 1 presents a preliminary 
analysis of the heating data obtained from Fire 1 and points out the need for an analysis 
of the body motions after the disturbance to assess their effect on the interpretation of 
the heating data. References 2 and 3 define the flight trajectory and atmospheric environ- 
ment and describe the mechanical hardware of Fire 1. 

Payload performance and stability instrumentation consisted of three rate gyros 
and three linear accelerometers. During reentry these instruments indicated small 
coning half-angles of about 1° which were initiated by the separation impulse. The anal- 
ysis of the body motions after the large disturbance at approximately 1666 sec was com- 
plicated by the fact that (1) the roll- and pitch- rate gyros became inoperative, (2) the 
linear accelerometers were not located on the center of gravity and their outputs were 
not continuous, but rather were commutated, (3) the telemetry experienced periodic noise 
disturbances, (4) the reentry package was lost by radar shortly after the onset of telem- 
etry blackout, and (5) the excursions of the yaw-rate-gyro-data oscillations exceeded 
calibration limits. 

There were seven intervals on the telemeter oscillograph records during reentry 
that were free of noise interference. The present study is limited to analysis of these 
seven periods. The primary purpose of this study is to analyze the motions and deter- 
mine the magnitude and frequencies of the angle of attack during the heating data period. 

In addition, the study was made to substantiate the trajectory parameters determined in 
the study of reference 2. 

In general, simulation techniques employing aerodynamic data obtained from wind- 
tunnel tests and solutions to the equations of motion as obtained from digital-computer 
programs were used to analyze the motion data. The ranges of velocity and altitude of 
the simulations were from about 11 430 meters per sec at approximately 63 400 meters 
to 2 800 meters per sec at 33 070 meters. 
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SYMBOLS 
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a x ,a y ,a z 
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C X 
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d H “ X o + AX 

d L = X Q - AX 
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F F F 
X' Y’ Z 


g 


normal acceleration, -A^, g units 

resultant acceleration in Y b -Z b body-axis plane, 
g units 

accelerations parallel to X b , Y b , and Z b body axes, 
respectively, g units 

resultant- or normal-force-curve slope 

resultant force coefficient in body-axis system 
force coefficient along X b body axis 

force coefficients due to asymmetry along Y b and Z b body axes, 
respectively 

resultant pitching- moment coefficient in Y b -Z b body-axis plane 
pitching- moment- curve slope 

moment coefficient about Y b body axis due to pitching velocity 

reference diameter, meters 
high damping rate, 1/sec 
low damping rate, 1/sec 
base of Napierian logarithm 

force components along X, Y fe , and Z b body axes, respectively, 
newtons 

acceleration due to gravity, meters/sec^ 
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h 


I 


I' = — , sec 2 
qSd 

i=V^ 


k 1 ; k 2 ,k 3 ,k 4 




m 

m' = sec 

qs 

q 

p,q,r 


Rj j R 2 

S 

t 

At n 

U,V,W 


u,v,w 


altitude, meters 

p 

mass moment of inertia about X b body axis, kilogram-meters* 

mass moment of inertia about Y b and body axes, 
kilogram- me ter s^ 


constants in angular- rate solution of linearized equations of motion 

moments about the and body axes, respectively, 

newton- meters 

mass, kilograms 

dynamic pressure, newtons/meter ^ 

angular velocities about the X^, Y^, and Z^ body axes, 

respectively, radians/sec 

initial magnitudes of rotating vectors of epicyclic motion 

reference area, meters^ 

elapsed flight time, sec 

time used in equations (1) and (2), t n - tj, sec 

inertial velocity components along X^, Y^, and Z^ body axes, 

respectively, meters/sec 

earth- relative velocity components along X b , Y b , and Z b body 
axes, respectively, meters/sec 
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V 

x b’ Y b’ z b 

x g’ Y g’ z g 

x 


x 


eg 


X 


o 


Vx> z x 


x y ,y y ,z y 


*z’\>\ 


y b ,Z b 

0 ! 

V 


r P 

r y 


earth- relative velocity, meters/sec 

orthogonal body-axis system with origin at center of gravity 

orthogonal gravity-axis system with origin at center of gravity 

longitudinal location of accelerometer measured from Fire 
station 0 (see fig. 2), meters 

longitudinal position of center of gravity measured from Fire 
station 0, meters 

longitudinal position of center of pressure of asymmetric force, 
meters 

displacements between longitudinal accelerometer and center 
of gravity parallel to X b , Y b> and Z b body axes, respectively, 
meters 

displacements between transverse accelerometer ( AyJ and center 
of gravity parallel to X b , Y b , and Z b body axes, respectively, 
meters 

displacements between normal accelerometer (A%J and center of 
gravity parallel to X b , Y b , and Z b body axes, respectively, 
meters 

coordinates along Y b and Z b body axes, respectively, meters 

angle of attack, deg 

angle of yaw, deg 

total angle of attack, deg 

earth- relative flight-path angle, deg 

earth- relative heading, deg 
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<P' 


AX 

v V v 2 

P 


W H 

W L 


body attitude angles relative to gravity-axis system, deg 
aerodynamic phase angle, deg 
nonrolling damping rate, per sec 
damping rate due to roll, per sec 

initial angular orientation of rotating vectors, radians 

air density, kilograms/meter^ 

high frequency, w 0 + Aco, radians/sec 

low frequency, w 0 - Aco, radians/sec 

natural pitch frequency, radians/sec 

basic oscillation frequency, radians/sec 

component of total pitch frequency resulting directly from roll, 
radians/sec 


= q + ir 
Subscripts: 

0 initial conditions (unless defined differently elsewhere) 

max maximum condition 

trim trim condition 

eg center-of-gravity value 

1 indicated instrument value 

A dot over a symbol denotes differentiation with respect to time. 
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TEST AND MEASUREMENTS 


In general, the Project Fire 1 vehicle consisted of an Atlas launch vehicle, a veloc- 
ity package containing an Antares II-A5 rocket motor, and the reentry package which con- 
tained the experimental apparatus. The reentry package was a blunt Apollo shaped body 
with a unique composite heat shield made up of two ejectable phenolic- asbestos layers 
sandwiched between three beryllium calorimeters which allowed three separate measure- 
ments of the heat pulse to be made during reentry. 

The Atlas put the Fire 1 spacecraft (combination of velocity and reentry packages) 
into a precise ballistic trajectory along the Atlantic Missile Range. After separation 
from the Atlas, the spacecraft was oriented to the proper Antares ignition attitude by the 
velocity-package control system and was spun up to approximately 18.85 radians/sec for 
spin stabilization. Sometime after burnout of the Antares rocket, the reentry package was 
cold- separated from the velocity package by a spring mechanism and reentered the earth's 
atmosphere at 11 574 meters per sec at an altitude of 121 920 meters and approximately 
8 028 kilometers downrange near Ascension Island. A more detailed description of the 
vehicle and flight events may be found in references 2 and 3. 

Flight Tests 

The sequence of significant flight events during reentry are listed in table 1. The 
reentry package separated from the velocity package at t = 1640.5 seconds. Approxi- 
mately 26 seconds later an unexpected disturbance occurred and large motions ensued. 

The trajectory was confirmed by radar until radar blackout at t = 1660.2 seconds. The 
latter part of the trajectory was extended by particle trajectory computer programs as in 
references 2 and 3. These extensions, however, were only used as a guide in the motion 
analysis. 

The motion instrumentation in the payload consisted of three rate gyros and three 
linear accelerometers. After the disturbance (t » 1666 sec), the amplitudes of the yaw 
rate exceeded the range of the yaw-rate gyro. The roll-rate gyro was inoperative and the 
pitch rate- gyro operated erratically. All the accelerometer data were commutated. The 
accelerometers also were located well off the center of gravity and consequently measured 
accelerations due to angular velocity and angular accelerations as well as the translatory 
accelerations of the vehicle. The ranges of all the instruments, together with the loca- 
tions of the accelerometers, are listed in table 2. Details of other reentry- package 
instrumentation such as calorimeters and thermocouples are given in references 1 to 3. 

A photograph showing the reentry portion of the telemeter oscillograph record (after the 
disturbance) is shown in figure 1. The seven data intervals free of noise interference are 
clearly shown in the figure. It may be noted that differences exist both in the amount of 
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noise and the time various events occur on the two playback records. The noise difference 
is probably due to the arbitrary nature of the telemetry decay. The time difference is 
slight, on the order of 0.03 second, and is probably a result of a time lag in the playback 
system. The data of the second playback are included in the analysis since it provided 
some data not available on the first playback. Since there was an overlapping of data in 
all cases where data from the second playback were used, the time difference was easily 
discerned and compensated for by shifting the time scale of the second playback data to 
make the overlapping portions congruent. 

Physical Characteristics and Ground Tests 

A sketch of the Project Fire 1 payload (reentry package) is shown in figure 2(a) 
and a sketch of a similar configuration which was tested in various facilities of the NASA 
wind-tunnel complex is shown in figure 2(b). In general, the overall linear and angular 
dimensions of the two bodies were similar; however, a prime difference was in the curva- 
ture and edge of the face of the wind-tunnel test model. The physical properties of the 
flight model are presented in tabular form in figure 2 and are illustrated graphically in 
figure 3. The wind-tunnel data employed for the reentry analysis are presented in 
figure 4 for angles of attack from 0 to 80°. These data are averaged from data obtained 
during tests at a Mach number of 4.65 (ref. 4) and from unpublished wind-tunnel data at 
Mach numbers of 5, 7.3, and 9. The moment center was determined by the equations 
shown in figure 2. The aerodynamic reference area and length were based on the model 
frontal area and diameter, respectively. 

Exploratory wind-tunnel tests were also made after the flight to determine whether 
reasonable amounts of uneven heat erosion of materials on the face of the reentry package 
and/or disturbances from unsymmetrical ejection of heat-shield segments could have 
induced aerodynamic forces and moments of the order to cause the aforementioned dis- 
turbance. Models were tested with various combinations of pie- shaped sections removed 
from their faces for this purpose. No significant change in the static aerodynamic coef- 
ficients was found. Only the wind-tunnel data used in the analysis of the motions of the 
reentry package are presented in this paper. 

METHOD OF ANALYSIS 

Simulation techniques were used in the motion analysis since incomplete rate- gyro 
data did not permit more direct methods of workup. A block diagram summarizing the 
techniques employed in the analysis of the Fire 1 data is shown in figure 5. In general, 
the simulations were made to match the on- scale data frequencies and amplitudes. In 
order to start the simulations, a separate study shown on the left side of the block dia- 
gram was first necessary to approximate required initial conditions of the variables not 
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measured. Once these conditions were determined, the simulation calculations were 
obtained with a six-degree-of -freedom computer program as shown on the right side of 
the block diagram. There were seven data periods. Each motion simulation was essen- 
tially independent of the other data periods. 

Motion Initial Conditions 

It should be noted in figure 4 that the pitching- moment coefficient is linear up to 
angle of attack of about 35°. Thus, a curve-fit program using the angular- rate-history 
solution to the linearized equations of motion was set up to approximate the pitch rate 
from measurements made by the yaw-rate gyro. The principal steps in the derivation of 
this solution are presented in appendix A. The derivation is based on that of reference 5, 
which presents the angle- of -attack- history solution for the linearized equations. By sub- 
stituting the angle-of-attack solution in the linearized force equation of reference 5 and 
solving for the angular rates, the following forms for the rate- gyro equations were 
obtained: 


r = Kj sin w^t + *2 cos co^t + Kg sin o>jjt + K 4 cos co H t (1) 

q = -K2 sin u> L t + Kj cos w L t + sin w H t - Kg cos co^t (2) 

where t is equal to At n . 

These equations are written in terms of the body-axis system for a body with 90° 
rotational and mass symmetry and negligible trim angles. The body-axis system is shown 
in figure 6. It is assumed that free- stream conditions are constant due to the brevity of 
each data period. All damping factors are neglected because of the high speed during the 
period of the analysis. Although the errors due to large angles of attack would have some 
effect on the accuracy of the linearized analysis, it was found that the utilization of the 
linear solution for obtaining approximate initial conditions for the exact calculations was 
justified on the basis that iterations could be performed if required in the final simula- 
tions. 

Basically, the linear solution is written in terms of the algebraic sum of two sine 
waves of different amplitudes, frequencies, and initial phase angles, and may be repre- 
sented by the epicyclic motion of two rotating vector arms revolving at different rates. 
There are six unknowns (kj, K 2 , Kg, K 4 , o> L , and to be determined in the anal- 
ysis of each data period. First, the frequencies u> L and o> H of the two rotating arms 
are computed from estimated values of relative velocity, dynamic pressure, roll rate, 
vehicle mass characteristics, and linear wind-tunnel aerodynamics (shown in eqs. (A9), 
(A10), (A14), (A22), and (A23)). These frequencies plus four values of yaw- rate data 
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chosen from within the first data cycle to be fitted were substituted in equation (1) to 
yield four simultaneous equations. The remaining unknowns or constants were evalu- 
ated from the solution of these four equations. The yaw rate is then generated from 
equation (1) and compared with a corresponding portion of the flight record. This 
procedure may be iterated by changing the estimated frequencies slightly, until a very 
close match is obtained. The time interval selected for curve fitting should be limited to 
keep the assumption of constant free- stream conditions valid. Finally, the constants were 
substituted in equation (2) to obtain the corresponding pitch-rate history. Equations (1) 
and (2) were differentiated to obtain the angular accelerations of the vehicle. With these 
angular-velocity and acceleration histories, all the remaining initial conditions, with 
exception of the trajectory quantities, necessary to start an exact solution using a six- 
degree- of -freedom computer program may be determined; a description of the procedure 
is presented in appendix A. 


Simulations 

The six- degree- of-freedom computer program is presented in reference 6; only 
those equations relative to generating the motion patterns are presented in appendix B. 
This program represents the exact solution to the equations of motion for six degrees of 
freedom. However, it is only exact to the degree required by the user. In the present 
instance, for example, the damping was neglected for the same reasons the damping was 
neglected in the curve-fit program. The exact program determines both translatory and 
rotational motions of a rigid body in three-dimensional space. Gravity effects and posi- 
tion were computed relative to an oblate spheroidal earth with a 1962 Standard 
Atmosphere (ref. 7). The variation with time of the center of gravity, mass, moments 
of inertia, reference areas and diameters, plus the complete nonlinear aerodynamic 
coefficients as presented in figures 2, 3, and 4, was programed. An option was added 
to this program to calculate the accelerations of an arbitrary point on the body. This 
option was added to simulate accelerometer data of instruments which had been located 
off the center of gravity of the reentry package (eqs. (B14), (B15), and (B16)). 

With given aerodynamic and mass characteristics, the high and low frequencies 
determined for each data interval by the curve-fit program were transformed to those 
quantities with which the exact program works in terms of — namely, q and p. This 
transformation was made through the use of equations (A24) and (A25). Subsequently, 
values of V and p were determined to satisfy the derived value of q. For simula- 
tion of the first data interval, V was obtained directly from the extension of the flight 
trajectory of reference 3 since it was believed that V was less sensitive to possible 
errors in that trajectory at the onset of reentry. The required air density p was deter- 

1 -2 

mined by solving a rewritten form of the usual q expression, - pV . By referring to 
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the U.S. Standard Atmosphere, 1962 (ref. 7), an initial altitude was then selected. A mini- 
mum of two simulation runs were made in the data period by using the same initial veloc- 
ity and altitude, but slightly different flight-path angles, and continued through the 
adjoining noise interference interval to the starting time of the next data period. This 
procedure was used to aid in determining the combination of velocity, altitude, and flight- 
path angle which would again satisfy the q requirement determined by the curve-fit 
program in the latter data period. This procedure was repeated for each of the 
remaining data periods. 

Once the data curves were matched, the computed histories of a, (3, rj, 9, and 
i\j were assumed to represent the flight motions and conditions and were used as a basis 
for a radiative -heating analysis of the Fire 1 payload (ref. 8). 

RESULTS AND DISCUSSION 
Basic Data 

Comparisons of the basic flight data with the numerical solution of the six-degree- 
of -freedom program for a portion of each data interval are presented in figures 7 to 13. 
The symbols in the rate-gyro curves represent selected points taken at constant intervals 
from continuous data, whereas the symbols on the acceleration curves indicate the com- 
mutated data. Flagged symbols represent data from a different playback of the record. 
The portion of the record presented represents virtually all the clear data in the interval. 
The comparisons include those of r, Ay p and plus one of q in the last period 

(t « 1687 sec) when it appeared the pitch- rate gyro might be operating properly. Also 
presented is a comparison of the generated yaw rate from the linear solution of the curve- 
fit program with flight and numerical computation values in the first period (see fig. 7(a)). 
This case is typical and more comparisons are not presented. A summary of the four 
data points selected from each data interval for the curve-fit program is presented in 
table 3 together with the resultant frequencies and amplitude constants obtained therefrom. 
Likewise, a summary of the initial conditions used in the exact solution of the computer 
program is presented in table 4. 


Simulations 

In general, the comparison in figures 7 to 13 shows that there is good agreement 
between the simulations and flight data as to amplitudes and frequencies. The simulation 
frequencies, in particular, are very close to those of the data. There is no direct way to 
compare the peaks of the amplitudes of either the rate- gyro or accelerometer data since 
the angular rates exceeded the calibration of rate gyros and the data of the accelerom- 
eters were commutated. Histories of the missing data beyond the limits of the yaw- rate 
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gyro were computed and are considered valid since close agreement was obtained with 
preceding and following on- scale data. The simulated values of Ay i and A z ^ pro- 
vided the only continuous accelerometer variation, without which it would have been most 
difficult to make valid fairings of the commutated flight data. The techniques used in the 
analysis of the accelerometer data were necessary for two reasons: (1) each accelerom- 

eter was located separately at different distances from the center of gravity and thus 
measured accelerations due to angular rates and accelerations (see appendix B) in addi- 
tion to the accelerations of the center of gravity and (2) since A z j and Ay^ were 
commutated, their recorded values necessarily occurred at different times; hence, direct 
computation of the resultant total acceleration Aj^ and the total normal-force coefficient 
C R would have been most difficult. 

The good agreement between the curve-fit-program curve and the numerical- 
solution curve (fig. 7(a)) is due primarily to the linearity of the C R and C m wind- 
tunnel data up to angles of attack of 35°. Trim forces and moments due to aerodynamic 
and mass asymmetries were neglected in the simulation of the accelerometer and rate- 
gyro data. The good agreement between the numerical solution and the flight data 
(fig. 7(a)), therefore, indicates that Vtrim the flight vehicle was small. 

Two significant parameters, Au> and co 0 , were isolated from the simulations and 
their variations with time are presented in figure 14. By neglecting the damping effects 
during small time intervals, it was possible to solve the expressions for these frequency 
parameters for roll rate p and the product C mo ,q. It was thought the roll rate would 
remain constant throughout the flight since the reentry package was a symmetrical body 
of revolution with no fins to provide rapid roll damping. However, it was found from the 
curve-fit program and validated by the exact program that p was initially 
18.85 radians/sec and generally decreased with time. This roll-rate deceleration was 
probably caused by such factors as skin-friction roll damping and the effects of possible 
heat erosion on the rolling moment. It may be noted that the dropoff of p with time was 
not uniform (see table 4) and was probably the result of sudden changes in mass charac- 
teristics and nose shape when the phenolic-asbestos heat shields were ejected 
(t ~ 1670 and 1676 sec). A study of the effects of inertia was not made since preliminary 
estimates showed that practical errors in 1^ and I would have only a small effect on 
the frequencies. Therefore, the aforementioned iterations were concentrated on finding 
values of rolling velocity and natural pitch frequency required to match the flight data. 

The variations Of pitching- moment- curve slope, based on the wind-tunnel data (see 
section entitled "Physical Characteristics and Ground Test") adjusted for reentry package 
center- of- gravity location referred to the wind-tunnel model, are shown in figure 15(a). 
The corresponding values of dynamic pressure employed in the simulation runs (table 4) 
are presented and compared with dynamic-pressure calculations from a ballistic 
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trajectory of reference 3 and a modified ballistic trajectory from reference 2 in 
figure 15(b). The modified trajectory was computed by using drag coefficients based on 
the average angle of attack determined from the present report. The differences shown 
in q from the three methods are not significantly large and may be accounted for by 
experimental errors in either the flight data or the wind-tunnel tests, or both. 

Angles of Attack 

The motion patterns in the inertial plane (0, \p) at selected flight intervals are pre- 
sented in figure 16. The representative altitudes and relative velocities for the starting 
time of each pattern are listed in table 4. Each interval is approximately 1 second long; 
hence, changes in altitude and velocity were small, and as a result, the pattern in each 
interval shown is generally symmetrical. However, the patterns are changing gradually 
with time. The patterns are typical of those for statically stable bodies which have roll- 

to-pitch inertia ratios greater than 1 > lj , inasmuch as such bodies always have below- 

resonance motions characterized by outside loops. The center of the inertial patterns is 
approximately the projected velocity vector defined by its two flight-path angles, y p 
and yy. The radial distance between the center of each pattern and the curves is approx- 
imately the resultant angle of attack. A comparison of the successive peaks in each pat- 
tern shows that the damping is small; the amount of damping or divergence indicated in 
successive patterns is due to the dynamic-pressure gradient since C m ^ = 0 in the 
calculations. 

The motion patterns in the relative plane (a, (3) at selected flight intervals are pre- 
sented in figure 17. These motions are symmetrical about the origin (a = 0, /3 = 0) since 
trim effects were zeroed. The radial distance from the origin to the curve is again 
approximately the total angle of attack. The effects of damping are similar to those noted 
previously. The difference between the relative and inertial patterns lies in the axis sys- 
tem to which each is referred. The relative motions are referred to the body-axis sys- 
tem which rotates with the body at its roll rate. The inertial motions are referred to the 
gravity-axis system which rotates only with respect to an inertial reference, as the body 
moves over the earth. 

The total angle-of-attack envelope is presented in figure 18 and shows the maximum 
and minimum bounds of its oscillations. Until t * 1666 seconds, the reentry package was 
coning with 1° half-angle cone thought to be caused by tipoff impulses at separation from 
the velocity package. A disturbance was experienced at t ~ 1666 seconds. (See fig. 1.) 
Based on the calculated coning-motion buildup thereafter, a resultant momentum impulse 
of about 19.4 newtons- meter- sec would have been required to increase the total angle-of- 
attack oscillations from 1° to about 33°. No attempt was made to analyze this disturbance 
further since the continuous data from the yaw-rate gyro experienced considerable noise 
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interference shortly after the disturbance. The upper bound of the 17 envelope decreased 
to a minimum of about 13° at approximately q max and then increased gradually as the 
dynamic pressure decreased. This angle- of- attack envelope is the one on which the modi- 
fied ballistic trajectory of reference 2 is based. 

CONCLUDING REMARKS 

An analysis of the flight motions of the Project Fire 1 payload was made by a 
digital- computer simulation of the rate- gyro and accelerometer data. This method of 
analysis was necessary since limited instrument data, consisting of continuous data from 
only the yaw-rate gyro, and commutated accelerometer data were available. The purpose 
of this analysis was to determine the angle- of-attack and motion histories of the reentry 
package during the experiment portion of the flight. The simulations were obtained from 
a six- degree- of-freedom computer program which numerically solved the equations of 
motion. The complete nonlinear wind-tunnel aerodynamics and calculated time histories 
of the mass characteristics of the body were used in the computer program. The reentry 
package experienced a disturbance at about 1666 seconds of elapsed flight time at an alti- 
tude of approximately 69 800 meters and a velocity of about 11 500 meters per sec. Just 
before the disturbance, the body was coning with a half-angle of about 1°. The results of 
the analysis indicate the following: 

(1) After the disturbance, the angle- of- attack envelope increased from 1° to a maxi- 
mum of about 33°, subsequently damped to a minimum of about 13° at maximum dynamic 
pressure, and increased somewhat thereafter as dynamic pressure decreased. While the 
factors which caused the disturbance could not be isolated, it was determined that a result- 
ant momentum impulse of about 19.4 newton- meters- sec would have been required to pro- 
duce this change in the magnitude of the motions. 

(2) The roll rate was initially 18.85 radians/sec and decreased slightly during the 
experiment portion of the flight. 

(3) The motion patterns of the reentry package were typical of those for bodies 
which are statically stable and have roll-to-pitch inertia ratios greater than 1. 

(4) Total angle -of -attack trim was small. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 17, 1965. 
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APPENDIX A 


ANGULAR-RATE SOLUTION OF THE LINEARIZED EQUATIONS OF MOTION 
AND APPLICATION FOR OBTAINING INITIAL CONDITIONS 
FOR SIMULATION OF FLIGHT 


Angular-Rate Solution of the Linearized Equations of Motion 

The assumptions made for linearizing the equations of motion for a spinning sym- 
metrical missile with respect to the body-axis system (fig. 6) are as follows: 

(1) Constant free- stream conditions exist. 

(2) The missile has constant rolling velocity. 

(3) Aerodynamic damping forces and gravitational forces are negligible in compari- 
son with forces due to angle of attack and angle of sideslip. 

(4) Small angle approximations are valid. 

(5) Magnus forces and moments are negligible. 

The equations of motion reduce to the following four expressions: 

" C N n ,^ S + C Y ® SF Y 
/3 + r - pa = mu = mu 


a - q + p/3 


-C^aqS + C^qS SFz 


mu 


mu 


(A2) 


ql - pr/l - I x ) = C m aqSd + C m M)qSd - C z qSx Q = 2M y (A3) 

\ / ot q\2V/ o 

i-I + pq(l - I x ) = -C m ^qSd + C m J|qS<l + C Y qSx a = 2M Z (A4) 

The solution of these linearized equations for total angle of attack t] is given in refer- 
ence 5 and by changing some symbols for convenience and clarity becomes: 


u>* 


V = ^trim + R l e e 


(d H +icu L ) 


t 


+ R 0 e 




(d L -iw H )t 


(A 5) 
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where in reference 5: 


r? = (3 + ia 

(A 6 ) 

d H = X 0 + AX 

(A 7) 

d L = X 0 - AX 

(A 8 ) 

GOr = Ct > 0 4- Ac 0 

(A9) 

co T = co - Aco 
L o 

(A10) 

l/ C N a C m q \ 
° 2 \m' r 7 

(All) 



AX 


1 

4a > 0 


fog 

I \m' 



(A13) 


Ao> 



(A14) 


Equations (Al) and (A2) may be combined by the use of complex numbers, and reduced, by 
letting = q + ir and m' = to the following form: 


• • . 0 ' C ^a n C Y 0 +iC Z 0 

77 + ip?y - + 


m' 


m' 


(A15) 


The angle-of-attack solution given in equation (A5) was substituted in the force equation 
(eq. (A15)) to yield the following solution for the angular rate: 


16 


APPENDIX A 


0 = q + ir = 


'N. 


p + co L -ildH + ^T 


a 


Rje 


d H t+i ( w L t+I 'l) 


'N 


p - co H - ild L + — 


a 


m’ 


R 2 e 




d T t 


C N a \ i ( C Y 0 + iC Z Q ) 


+ \ P " 1 l^ftrim 


m T 


(A16) 


By using the relation e ir = cos t + i sin t, the time histories of the pitching rate q and 
yawing rate r can be obtained from equation (A16) as follows: 


q + ir = R-i e 


d R t 


N r 


(P + "l) 008 ^ + ^l) + \ d H + -^/ sin ( w L t + v l) 


+ I^e 




(p - “h) 008 ^ - v 2) - ( d L + " v *) 


/ '"'Z 0 \ . 

+ (P^trim + m i “trim “ m i J + 1 


R l e 


d H t 


(p + w L )sin(o> L t + vij 


d H + 


'N. 


_^)cos(a, L t + Vl ) 


+ R2e 


d L f 


"(P “ “h) 810 ^ • v 2) 


'N 


d L + -^j C0S (^ - v l) 


( U N \ 

+ (p^trim " m' ^trim + m ' ) 


(A17) 


By using fundamental trignometric identities, algebraically collecting terms, and sepa- 
rating the real from the imaginary parts, the following expressions for q and r may 
be written: 


r = e 


d H t 


Ri 


N„ 


(p + u> L )cos V X + \d H + ——~jsin V X 


sin Wj_t 


d H tl 

+ e Ri 


(p + cu L )sin v x - (d R + -^yjcos 


cos a> L t 


(Equation continued on next page) 
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+ e 


d h t 


Rr 


N fl 


-(P - w h) cos v 2 - ( d L + -^r/ sin v 2 


'sin Wjjt 


+ e 


d Lt 


Rr 


(P “ w H) sin v 2 ~ ( d L + -J^) COS v 2 


COS COjjt 


/ C N a C Y 0 \ 

+ (P^trim " nV ^trim + m' / 


(A18) 


drrt 

q = e H (Rj 


-P 


/ C N«\ 

+ cousin ^ + \d H + —Jcos ^ 


sin co^t 


+ eHt, 


-N, 


(P + W L) C0S V 1 + \ d H + -^T/ Sin v l 


cos u>j_t 


d T t\ 

+ e L R r 


/ C N c \ 

- cousin v 2 - \d h + —Jc os v 2 


i sin u>jjt 


+ e dL,t |R 


N, 


(P - w H )cos V2 + [d L + —£Lj sin v 2 


[cos a»jjt 


N a __ C Z 0 \ 


+ V^trim + m f ^rim “ m T 


m T 


(A19) 


where the last terms in parentheses in equations (A18) and (A19) are the r^ r j m and 
qtrim terms, respectively. By simplifying the constant terms, the equations reduce to: 

r - r t r j m = e 11 (Kj sin a> L t + Kg cos o> L tj + e^^^Kg sin Wjjt + cos Wjjtj (A20) 
q - Qtrim = e H (-Kg sin + K 1 cos w i_t) + e L ( K 4 s t n W H* " ^3 cos “H*) (A21) 
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Equations (A20) and (A21) are the equations which define the quantities r and q meas- 
ured by the rate gyros. For the present application the time interval selected for each 
analysis was less than 1 second of flight time; thus, for these short time intervals 


Also, examination of the basic data did not clearly indicate any significant effect, if any, 
of trim on the rates, and therefore, the q|. rim and r trim terms were neglected. 
Equations (A20) and (A21) then reduce to the form used in the curve-fit program (eqs. (1) 
and (2)), namely: 


r = Kj sin w L t + K 2 cos o> L t + Kg sin co H t + K 4 cos co H t 


q = -K 2 sin co L t + Kj cos co L t + K 4 sin co H t - Kg cos u> H t 
Also, equation (A12) may be reduced to: 


where 


and 



(A22) 



“C mQ qSd 

i 


(A23) 



Simulation Initial Conditions 

Given q Q , r Q , c^, r Q , o> L , and at an arbitrary starting time within the 

data interval from the curve-fit program, the remaining initial conditions necessary to 
start an exact solution using a six-degree-of-freedom computer program to be deter- 
mined are: the dynamic pressure q Q , the roll rate p Q , the aerodynamic phase angle 
cp' Q , the total angle of attack ri 0 , and the body attitude angles 9 Q> \p Q , and cp Q . The 
initial dynamic pressure is determined by combining equations (A9), (A10), (A22), and 
(A23) to yield 
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-C m Sd 
m a 


+ w L 


<■(•4 


(A24) 


where the roll rate p Q is determined by solving equations (A9), (A10), and (A14) to give: 


W H ~ W L 


(A25) 


When equations (A3) and (A4) are modified to permit the use of nonlinear aerodynamics, 
the initial aerodynamic phase angle (p' Q , the orientation angle of the plane of the velocity 
vector with respect to the body- axis system, may be determined in the following manner: 


rl + pq^I - I x j = 2M Z = -C m sin cp' qSd 


(A26) 


ql - pr^I - I x ^ = 2 My = C m cos <p' qSd 


(A27) 


Dividing equation (A26) by equation (A27) yields 


-l|-|fl pq(! - I x ) 


<P'o = tan 4 -L 2U. 

ql - pr(l - I x ) 


(A28) 


Since C mQ) is essentially constant up to an angle of attack of approximately 35°, the 
initial total angle of attack may now be obtained from 


_ C m,o 

''O c 

'-'tn 


(A29) 


where 


m ’° \-sin <p' qSd/ Q \cos <p' qSd/ Q 


Finally, the body attitude referred to the gravity- axis system (fig. 6) is determined with 
the assumption that the initial roll position <p Q may be arbitrarily selected since it has 
little effect, if any, on the relative motion patterns of the body and affects only the orien- 
tation of the inertial motion patterns. The body pitch and yaw angles 0 and ip, respec- 
tively, may now be obtained by trigonometrically adding Vo anc * <Po to the pitch and 
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yaw flight- path angles y and y respectively, from the trajectory to yield: 

r J 


_ Jsin yp cos rj + sin rj cos c pQ -JcosPr] + sin Prj cos %<pQ - sin^yp^ 


0 Q = sin' 


cos^r] + sin^rj cos^ 


and 


where 


*o=r y ,o- sin ' 


Wsin r\ sin (Pq\ 

\ C0S V )o 


V G,o “ (”' - ”)o 


(A31) 


(A3 2) 


(A33) 
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EQUATIONS OF MOTION USED IN THE SIX-DEGREE-OF-FREEDOM 

COMPUTER PROGRAM 

The equations used in the simulations of this analysis may be found in reference 9, 
with the exception that the trajectory calculations were made with respect to an oblate 
spheroidal earth rather than the spherical earth of reference 9. These equations are a 
part of a general computer program which is described in detail in reference 6. Only the 
equations directly related to the simulation of rate and accelerometer data and to the 
generation of corresponding relative and inertial motion patterns are presented herein. 

To simulate the accelerometer data of instruments located off the center of gravity, the 
accelerometer transfer equations from reference 10 were added to those of the general 
computer program. The motion equations of the program are written for a rigid body 
with variable mass and determine both the translatory and rotational motions of the center 
of gravity of the body in three-dimensional space. The Standard Atmosphere, 1962 
(ref. 7), was programed. Unlike the curve-fit program of appendix A, there were no 
restrictions on free- stream conditions, gravitational forces were considered, and allow- 
ance was made for large angles and nonlinear aerodynamics. Roll rate, however, was 
held constant within each interval, although it was permitted to change from one data 
period to another. The forces and moments due to damping, Magnus effects, and aero- 
dynamic and mass asymmetries were neglected. With the preceding conditions, the 


equations were reduced to the following: 

m(u + wq - vr) + mg sin 8 = mgAx >C g = C^qS (Bl) 

m(v + ur - wp) - mg cos 8 sin <p = mgAy C g = -C R sin <p' qS (B2) 

m(w + vp - uq) - mg cos 8 cos <p = mgA Z)Cg = -C R cos <p' qS (B3) 

pl x = 0 (B4) 

ql - pr^I - I x ^ = C m cos <p' qSd (B5) 

rl + pq/l - I x \ = -C m sin cp' qSd (B6) 
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The earth- relative velocity V, was calculated by vectorially subtracting the rota- 
tional velocity of the earth from the inertial velocity of the body. The angles of attack 
and yaw of the vehicle were computed in terms of V and its components along the body 
axis, u, v, and w, as follows: 


a = tan"*- ^ 

(B7) 

u 


|3 = sin -1 

(B8) 


V 

The total angle of attack and its phase angle may be defined in terms of the components 
angles, a and /3, as follows 


t] = cos~l(cos a cos /3) 


(B9) 


<p' = tan" 1 iSHLS. (BIO) 

sin a 

The trajectory equations of the program are the usual three-dimensional expressions and 
are not repeated here. 

The attitude of the vehicle with respect to the gravity-axis system is defined by 
Euler angles (see fig. 6) and are calculated from the following equations: 


0 =0 O + J 

(q cos <p - r sin <p - K^idt 

(Bll) 

1) 

o 

+ 

'/r cos <p * q sin <p + -Y* 
\ cos 6 7 

(B12) 

<P = <p O 

+ C(p + i/' sin 6 + K<p)dt 

(B13) 


where Kg, K^, and are the body-axis components of the angular velocity of the 

local gravity-axis system moving with the vehicle with respect to an inertial reference. 
(See ref. 9.) 

Finally, the equations for determining the off-center-of- gravity accelerations 
(ref. 10) were obtained from: 

A x, i = A x, eg - 1^ 2 + r2 ) + y# - p<j) “ z x(4 + r P)] < B14 ) 
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A Y, i = A y,cg - |(- x y (f + pq) + y y ( r2 + p2 ) + z y (p " qr) l 

A Z, i = A Z,cg - |[*zW - rp) - y z (P + qr) + z z (p 2 + q 2 )] 


(B15) 

(B16) 
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TABLE 1.- SEQUENCE OF EVENTS DURING REENTRY 


Event 

Elapsed flight time, t, sec 

Reentry-package separation 

1640.5 

Reentry commences (h = 121 920 m) 

1647.4 

Begin telemeter blackout 

1653.9 

Begin C-band radar blackout 

1660.2 

Disturbance 

1665.94 

First heat- shield ejection (signal) 

1669.6 

Second heat- shield ejection (signal) 

1676.6 

End telemeter blackout 

1686.8 

Reentry-package impact 

1965.7 
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TABLE 2.- MOTION-INSTRUMENTATION RANGES 


PLUS ACCELEROMETER LOCATIONS 


Motion instrument 

Range 

Location 

x, m 

(Fire sta.) 

y b > m 


Axial accelerometer, A x 

0 to -120g 

0.4550 

0.0475 

-0.0393 

Transverse accelerometer, Ay 

+6g to -6g 

.5105 

0 

-.0315 

Normal accelerometer, A z 

+6g to -6g 

.5105 

.0546 

0 

Roll- rate gyro, p 

0 to 35 radians/sec 




Pitch- rate gyro, q 

±3 radians/sec 




Yaw- rate gyro, r 

±3 radians/sec 
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TABLE 3.- SUMMARY OF INPUT AND RESULTS OF CURVE-FIT PROGRAM 

FOR SEVEN DATA INTERVALS 


Date period . 

Time interv 
sec . . . . 

ti . . . 


I 

n 

in 

IV 

V 

VI 

VII 

1667.64 

1670.36 

1673.06 

1676.21 

1679.61 

1683.17 

1686.85 

1668.38 

1671.08 

1673.84 

1676.95 

1680.38 

1683.91 

1687.68 

1667.64 

1670.36 

1673.06 

1676.21 

1679.61 

1683.17 

1686.85 

1667.70 

1670.40 

1673.08 

1676.23 

1679.63 

1683.22 

1686.92 

1667.92 

1670.62 

1673.26 

1676.30 

1679.67 

1683.26 

1686.99 

1668.05 

1670.65 

1673.28 

1676.32 

1679.71 

1683.32 

1687.05 


-2.793 -2.863 


2.915 


2.498 


0 

.02 

.20 

.22 

1.116 

-2.175 


0 

.02 

.09 

.11 

2.874 

-2.323 


0 

.02 

.06 

.10 

2.862 

-2.072 


0 

.05 

.09 

.15 

-2.974 

.617 


r 3 

2.887 

1.991 

-1.874 

-1.911 

-2.846 

-2.784 

-2.69 

r 4 

-2.649 

-2.793 

2.666 

2.482 

.830 

2.478 

2.98 

u, L .... 

9.6 

14.475 

22.100 

27.700 

32.70 

27.85 

19.03 

"!!•••• 

23.74 

26.055 

35.36 

43.300 

45.16 

40.06 

30.72 

Ki .... 

6.767 

8.318 

-4.877 

-3.385 

-.704 

-4.645 

-2.02 

k 2 .... 

-1.279 

-5.167 

5.486 

.625 

5.129 

-1.947 

-2.84 

k 3 .... 

.528 

1.285 

-2.656 

-3.325 

-5.490 

5.620 

4.61 

K4 .... 

-1.514 

2.304 

-4.370 

2.252 

-2.267 

-1.027 

.14 


At n - t n - tj. 

tr represents yaw rate data points used in eq. (1). 




-5.3i 
















General dimensions and mass characteristics. 

Figure 2.- Physical properties of flight and wind-tunnel models. Linear dimensions are in meters. 








'e 5.- Block diagram summarizing simulation techniques used in analysis of Fire 1 flight 
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Flight data 
Numerical solution 
Curve-fit program 



(a) Yaw rates. 

Figure 7.- Comparison of computed data with corresponding measured flight data. 1667.64 to 1668.38 seconds. Crossed symbols indicate data points 

selected for curve-fit program. 
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1671*0 
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t, seconds 

(b) I ndicated accelerations. 


Figure 8.- Concluded. 



radians/ second 



(a) Yaw rates. 

Figure 9.- Comparison of computed data with corresponding measured flight data. 1673.06 to 1673.84 seconds. Crossed symbols indicate data points 

selected for curve-fit program. 
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(b) Indicated accelerations. 
Figure 9.- Concluded. 
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radians/ second 



I 676.2 .3 A .5 .6 .7 .8 .9 1677 . C 


tj seconds 


(a) Yaw rates. 

Figure 10.- Comparison of computed data with corresponding measured flight data. 1676.21 to 1676.95 seconds. Crossed symbols indicate data points 
selected for curve-fit program. Flagged symbols represent data from a different playback. 
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radians/ second 


O Flight data 



(a) Yaw rates. 

Figure 11.- Comparison of computed data with corresponding measured flight data. 1679.61 to 1680.37 seconds. Crossed symbols represent data points 
selected for curve-fit program. Flagged symbols represent data from a different playback. 
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t, seconds 


(a) Yaw rates. 

Figure 13.~ Comparison of computed data with corresponding measured flight data. 1686.85 to 1687.68 seconds. Crossed symbols indicate data points 

selected for curve-fit program. 
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(c) Pitch rates. 


Figure 13.- Concluded. 
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Figure 17.- Comparison of angular-motion patterns in relative reference plane. 
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